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ABSTRACT 

We have derived azimuthally-averaged radial iron abundance profiles of the X-ray gas contained within 
12 clusters of galaxies with redshift 0.03 < z < 0.2 observed with BeppoSAX. We find evidence for a 
negative metal abundance gradient in most of the clusters, particularly significant in clusters that possess 
cooling flows. The composite profile from the 12 clusters resembles that of cluster simulations of Metzler 
& Evrard (1997). This abundance gradient could be the result of the spatial distribution of gas-losing 
galaxies within the cluster being more centrally condensed than the primordial hot gas. Both inside and 
outside the core region, we find a higher abundance in cooling flow clusters than in non-cooling flow 
clusters. Outside of the cooling region this difference cannot be the result of more efficient sputtering 
of metals into the gaseous phase in cooling flow clusters, but might be the result of the mixing of low 
metallicity gas from the outer regions of the cluster during a merger. 

Subject headings: cooling flows — galaxies: clusters — intergalactic medium — X-rays: galaxies 



1. INTRODUCTION 

The presence of metals in the hot gas contained within 
clusters of galaxies provides important clues as to how met- 
als are deposited into the intracluster medium (ICM) from 
the constituent galaxies of the cluster. The two most likely 
mechanisms for injecting metals into the ICM are through 
ram pressure stripping of the gas from galaxies by the ICM 
(e.g., Gunn & Gott 1972; Gaetz, Salpeter, & Shaviv 1987) 
and galactic winds (e.g., Ostriker & Yahil 1973; Metzler 
& Evrard 1994). Recent work has suggested that negative 
radial abundance gradients are a common feature in clus- 
ters. Analysis of ASCA data have indicated that in some 
clusters the metal abundance declines from ~40%-60% of 
solar in the center to ~20% of solar at a distance of 0.5-1 
Mpc (Ikebe et al. 1997; Ezawa et al. 1997; Sarazin, Wise, 
& Markevitch 1998; Finoguenov, David, & Ponman 1999; 
Dupke & White 2000a,b) 

Most of the above studies, however, dealt with groups or 
cooler clusters. The large, energy-dependent Point Spread 
Function (PSF) of ASCA preferentially scatters hard pho- 
tons, making spatially-resolved spectroscopy problematic 
for hot clusters. Whereas the problem is less severe for 
cooler (< 5 keV) clusters, substantial errors in the radial 
temperature and abundance profiles can occur if the PSF 
is not dealt with properly (see Takahashi et al. 1995). As 
such, little observational work has been done on the radial 
abundance profiles of hotter clusters. There has also been 
correspondingly little theoretical work performed on the 
predicted abundance profiles of clusters from hydrodynam- 
ical simulations. White (1999) has derived PSF-corrected 
abundance profiles for a large sample of clusters observed 
with ASCA, but the deconvolution technique used to cor- 
rect for the PSF unfortunately led to rather large uncer- 
tainties in the derived abundance profiles. As a result, 
only weak constraints could be placed on the presence of 
abundance gradients from these data. 

'^ Chandra Fellow. 



BeppoSAX is better-suited to perform spatially-resolved 
spectroscopy on hot clusters. The PSF of BeppoSAX is 
one-half that of the ASCA CIS, and more importantly, is 
only weakly dependent on energy. The abundance pro- 
files of five higher temperature clusters have already been 
derived from BeppoSAX, finding at least marginal evi- 
dence for abundance gradients in A2029, A2256, A3266, 
and PKS0745-191, with no evidence for a gradient in 
A2319 (Molendi & De Grandi 1999; Molendi, De Grandi, 
& Fusco-Femiano 2000; De Grandi & Molendi 1999a,b; 
Molendi et al. 1999). In this Paper, we increase this sam- 
ple size by analyzing BeppoSAX archival data for 12 clus- 
ters of galaxies, and derive radial abundance profiles for 
each cluster, nine of which have temperatures greater than 
5 keV. Throughout this paper, we assume Hq = 50 km s~^ 
Mpc~i and go = 0.5. 

2. SAMPLE AND DATA REDUCTION 

From the i?e ppo6'^V Science Data Center (SDC) ar chive 
(available at [http://www.sdc.asi.it/sax_niain.htmlD we 
have obtained data for 12 clusters of galaxies. Within this 
sample, eight of them possess cooling flows (A85, A496, 
A1795, A2029, A2142, A2199, A3562, and 2A0335+096) 
and four do not (A2163, A2256, A2319, and A3266). We 
define a cluster as having a cooling flow if the cooling rate 
is greater than 20 Mq yr~^; see Peres et al. 1998 and 
White, Jones, & Forman 1997. All of the clusters are at 
low redshift (0.03 < z < 0.09) except for A2163, which 
has a redshift of z = 0.203. We analyze data taken with 
the Medium Energy Concentrator Spectrometer (MEGS) 
onboard BeppoSAX. The MEGS (see Boella et al. 1999 
for details) is composed of two identical gas scintillation 
proportional counters (three detectors before 1997 May 9) 
that are sensitive in the 1.3-10.5 keV energy range. The 
event files for all 12 clusters were subjected to the standard 
screening criteria of the BeppoSAX SDC. 



Since we are interested in deriving the radial abundance 
profiles for these clusters, it is important to account for 
scattering owing to the (PSF) of the MECS instrument. 
Fortunately, the detector + telescope PSF of BeppoSAX 
is nearly independent of energy. This is because the Gaus- 
sian PSF of the MECS detector improves with increas- 
ing energy, while the PSF of the grazing incidence Mirror 
Unit degrades with increasing energy (D'Acri, De Grandi, 
& Molendi 1998), leading to a partial cancellation when 
these two effects are combined. 

To correct for the PSF we have used the routine effarea, 
available as part of the SAXDAS 2.0 suite of BeppoSAX 
data reduction programs. The program ejfarea convolves 
the surface brightness profile of the cluster (which has been 
determined using ROSAT PSPC data; see Mohr, Math- 
iesen, & Evrard 1999 and Ettori & Fabian 1999) with the 
PSF of the MECS to determine the amount of contamina- 
tion in the spatial region in question at each energy from 
other regions of the cluster. Thus, an energy-dependent 
correction vector is formed. When this correction vector 
is multiplied by the observed spectrum, the PSF-corrected 
energy spectrum is obtained for the region. In practice, 
this information is incorporated into the auxiliary response 
file (the .ar/ file), which is subsequently used in the spec- 
tral fitting. This task also corrects for vignetting. A more 
complete description of the task is given in Molendi (1998) 
and D'Acri et al. (1998). This correction does not appear 
to affect the spectrum greatly; D'Acri et al. (1998) and 
Kaastra, Bleeker, & Mewe (1998) found only small changes 
between the uncorrected and corrected temperature pro- 
files for Virgo and A2199, respectively. In addition, D'Acri 
et al. (1998) found that the correction vectors amounted 
to 5% or less for energies above 3 keV outside of the inner- 
most bin (the innermost bin loses some flux via scattering 
but gains very little from photons scattered in from greater 
off-axis radii). Our inner three spatial bins are identical 
in angular extent to the bins of D'Acri et al. (1998), so we 
assume that our correction vectors are as small as theirs. 
In addition, we find only a small difference between our 
corrected and uncorrected abundance profiles for our sam- 
ple. 

We extracted spectra in four annular regions for each 
cluster, with inner and outer radii of 0' — 2', 2' — 4', 4' — 6', 
and 6' — 9'. At 9' the telescope entrance window support 
structure (the strongback) strongly absorbs X-rays, mak- 
ing spectroscopy difficult in this area of the detector. In 
addition, for off-axis angles greater than 10', the depar- 
ture of the PSF from radial symmetry becomes noticeable 
(Boella et al. 1997). With this in mind, we ended the radial 
profiles at 9'. At this radius, the profiles extended to 17%- 
33% of the virial radius, where Vmriai = 3.9 (T/10 keV)^/^ 
Mpc (see Evrard, Metzler, & Navarro 1996), for all clus- 
ters except A2163, for which the profile extended to 55% of 
T virial ■ Wc also extracted one global spectrum (0' — 9') and 
also two spectra covering the two regions — 0.075ri,iriai 
and 0.075 — QATHrmriai- Background was obtained from 
the deep blank sky data provided by the SDC. We used 
the same region filter to extract the background as we did 
the data, so that both background and data were affected 
by the detector response in the same manner. The energy 
channels were regrouped to contain at least 25 counts. 

For each cluster, XSPEC Version 11.0 was used to fit the 
spectrum. The MECS2 and MECS3 (and MECSl when 



available) data were fit separately, but with the same tem- 
perature, metallicity, and normalization. We assumed a 
MEKAL model with an absorption component fixed at the 
Galactic value, allowing the temperature and metallicity 
to vary. For the iron-to-hydrogen ratio we use a value of 
4.68 X 10~^ (Anders & Grevesse 1989) since this is the value 
assumed for most previous abundance determinations, al- 
though it should be noted that more recent determinations 
of the Fe/H ratio point to a value of 3.24 x 10"^ (Ishimaru 
& Arimoto 1997). Because of a systematic shift of 45-50 
eV at 6.6 keV in the MECS channel-to-energy conversion, 
we have also allowed the redshift to vary (F. Fiore, private 
communication). This systematic shift was evident in our 
sample; when the redshift was allowed to vary, the mea- 
sured redshift was less than the optically-determined red- 
shift in all 12 clusters, and inconsistent with the optically- 
determined redshift at the 90% confidence level for eight of 
them. A modest decrease in the reduced x^ also occurred 
for most of the clusters when the redshift was allowed to 
vary. However, freeing the redshift did not affect the val- 
ues obtained for the temperature and metallicity by more 
than 5% (and in most instances much less), so the results 
do not depend significantly on this shift. We also fit only 
the 3.0-10.5 keV part of the spectrum because of excess 
flux found below 3.0 keV that might be the result of uncer- 
tainties in the calibration of the MECS instruments; see 
Irwin & Bregman (2000) for a more complete description. 
All quoted errors are 90% confidence levels for one inter- 
esting parameter (Ax^ = 2.71) unless otherwise noted. 
Below we present the results from the abundance analysis. 
The results of the temperature analysis are given in Irwin 
& Bregman (2000). One additional cluster (A3562) not 
analyzed in Irwin & Bregman (2000) is presented here. 

3. ABUNDANCE OF THE HOT GAS 

3.1. Global Abundances 

For hot clusters it is the Fc-Ka line that primarily de- 
termines the metal abundance measurement. The fits to 
the global spectra were adequate, with xt < 1-20 in all 
cases. For the eight clusters harboring cooling flows, the 
weighted average emission- weighted global abundance was 
0.37±0.04 (the error is given by the weighted standard de- 
viation for the subsample). The four clusters lacking cool- 
ing flows had an average emission-weighted global abun- 
dance of 0.27 ± 0.02. The ratio of cooling flow to non- 
cooling flow global abundance is 1.4 ± 0.2, in good agree- 
ment with the value of 1.6 obtained by Allen & Fabian 
(1998) with a sample of 30 clusters observed with ASCA 
when they used a single component thermal model. 

3.2. Radial Abundance Profiles 

The radial metal abundance profile for each cluster is 
shown in Figure |^. Four of the clusters have had abun- 
dance gradients determined with ASCA, and are generally 
consistent within the errors with the profiles determined 
here: A496 (Dupke & White 2000b; Finoguenov et al. 
1999), A2199 (Finoguenov et al. 1999), A2029 (Sarazin 
et al. 1998; Finoguenov et al. 1999) and 2A0335-I-096 
(Kikuchi et al. 1999). In addition, analysis of this same 
BeppoSAX data for A2029, A2256, A2319, and A3266 by 
other authors are in agreement with our results (Molendi & 
De Grandi 1999; Molendi et al. 1999; De Grandi & Molendi 
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Fig. 1. — Radial abundance profiles for the 12 clusters in the sample with 90% confidence levels, derived from spectral fitting in the 3.0-10.5 
keV energy range. Solid lines represent abundances corrected for the PSF of BeppoSAX and the dotted lines are uncorrected for the PSF, 
except for A2256 where the dotted line represents the BeppoSAX profile derived by Molendi et al. (2000) from a longer observation (their 
error bars have been multiplied by 1.65 to convert 68% confidence levels to 90% confidence levels). 



1999), with the exception of A2256 for which Molendi et 
ah (2000) found a significantly decreasing profile whereas 
ours is flat. However, our inability to find an abundance 
gradient is the result of the short observing time of the 
data we used; the total observing time of the Molendi et 
al. (2000) observation was 2.6 times longer than ours, and 
they were able to place much tighter constraints on the 
abundance than we were. In both studies, the la error 
bars overlapped one another. 

In general, the abundance profiles decline with radius. 
Such a trend was also found in clusters observed by ASCA 
(Finoguenov et al. 1999; Kikuchi et al. 1999; Dupke & 
White 2000a), although there are notable exceptions to 
this trend such as A1060 (Tamura et al. 1996). These pre- 
vious studies focused primarily on cooler clusters {kT < 5 
keV). This study confirms the ubiquity of negative abun- 
dance gradients in hotter clusters (nine of our 12 clusters 
have kT > 5 keV), as illustrated in Figure |[ where all 12 
abundance profiles, normalized to the global abundance 
for each cluster, are plotted versus radius in units of the 
virial radius. A strong negative gradient in the abundance 
profile is evident out to 25% of the virial radius. However, 
there are a few points quite discordant with this trend at 
large radii. The clusters responsible for these discordant 
points are A1795, A2142, and A2163. 

We investigated the possibility that the use of single 
component temperature models in the central regions of 



cooling flow clusters has biased the determination of the 
abundance in these regions. The addition of a cooling flow 
component did not lead to a significant improvement in 
the fits, and the cooling rate was only weakly constrained. 
This is a result of excluding data below 3.0 keV, where the 
effect of cooling flows on the spectra is most pronounced. 
When we fixed the cooling rate at values determined by 
other satellites, the best-fit abundance value was actually 
slightly higher than in the single component case by a few 
percent, indicating that the central enhancement in the 
abundance is real. 

We have fit a linear model of the form {Z/Z„iean) = 
a — h(r JTmriai) to the data. While including all values, 
the fit is not very good, with x^ — 72.2 for 47 degrees of 
freedom, and can be excluded at > 98% confidence level. 
Removal of the four discordant data points yielded an ex- 
cellent fit, with x^ = 46.1 for 43 degrees of freedom for 
values of a = 1.212 and h = 2.705 (2.149-3.250 at 90% 
confidence) . Although it is not clear why four of the clus- 
ters show a ~2cr departure from the observed trend outside 
~25% of the virial radius, the derived slope is nonetheless 
an accurate representation of the data inside of 25% of 
the virial radius. Excluding non-cooling flow clusters in 
addition to the discordant points led to the same best-fit 
slope. 

3.3. Comparison To Theoretical Models 
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Fig. 1 - continued. 
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Surprisingly little theoretical work has been performed 
to date to quantify the behavior of abundance gradients 
in clusters. One such study is that of Metzler & Evrard 
(1997), who used gas dynamical simulations which in- 
cluded the effects of galactic winds to generate an ensem- 
ble of 18 realizations of clusters spanning a wide range of 
temperatures. Assuming that the iron present in the gas 
resulted from the ejection of metals from galaxies, a com- 
posite iron abundance profile for the 18 realizations was 
determined. It was found to decrease significantly with 
increasing radius. This resulted from the spatial distribu- 
tion of the galaxies being more centrally condensed than 
the original primordial X-ray gas. Since the ejected gas 
traces the galaxies, an abundance gradient is created rel- 
ative to the primordial gas. 

The composite iron abundance profile from the 18 three- 
dimensional realizations of Metzler & Evrard (1997) is 
shown superposed on the BeppoSAX data in Figure as 
a solid line. The composite profile was normalized as- 
suming an average global abundance of 40% of the solar 
value. Metzler & Evrard (1997) scaled each cluster in ra- 
dial units of ri7o before merging the profiles, where rno 
is the radius within which the mean density is 170 times 
the critical density. This is very close to our adopted value 
of Tmriai ■ The normalized abundance profile of Metzler & 
Evrard (1997) matches the data very well, ranging from a 
value of 1.3 in the center to about 0.6 at 25% of the virial 
radius. 

We note that recent work by Finoguenov et al. (1999) 
and Dupke & White (2000b) has shown that significant 
radial gradients exist in the elemental abundance ratios in 



clusters. This would argue against the idea that a differ- 
ence in the spatial distribution between the galaxies and 
the gas is responsible for the observed abundance gradi- 
ent if it is assumed that both Type la and Type II el- 
ements were released into the ICM under the same con- 
ditions. The observational evidence implies that Type la 
and Type II elemental contamination is not the same in the 
inner and outer parts, with Type la ejecta being dominant 
in the center. Dupke & White (2000b) propose that the 
observed abundance ratios can be explained by a phase 
of early, vigorous protogalactic wind driven by Type II 
supernovae that spread metals uniformly throughout the 
cluster, followed by a phase of less vigorous Type la super- 
novae wind that was suppressed by the central dominant 
galaxy in the center of the cluster. Thus, the measured 
iron abundances originate from a combination of Type II 
supernovae (having no radial gradient) and Type la super- 
novae (having a radial gradient). The gradient we measure 
is the linear combination of these two effects. 

3.4. Cooling Versus Non-cooling Flow Clusters 

Inspection of the abundance profiles of the eight cool- 
ing flow clusters reveals a radially decreasing trend in all 
of them. The gradient is quite significant in A85, A496, 
A2199, A3562, and 2A0335-t-096, and less significant in 
A1795, A2029, and A2142. This seems to be a universal 
trend in cooling flow clusters since all nine cooling flow 
clusters in the sample of Finoguenov et al. (1999), as well 
as all three clusters in the samples of Kikuchi et al. (1999) 
and Dupke & White (2000a) show the same trend, at least 
marginally. No clear trend exists among the non-cooling 
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flow clusters. The profile of A2256 is flat, although a de- 
cline with radius cannot be ruled out owing to the large 
errors (a decline was seen with a longer BeppoSAX ob- 
servation by Molendi et al. 2000; see Figure ||). A2319, 
A3266, and A2163 show some evidence for a negative gra- 
dient although the profiles are formally consistent with a 
flat profile. Results with ASCA for clusters with little or 
no cooling flows are also mixed. Whereas Coma, A401, 
A1060, and A2670 show a flat abundance proflle, A399 
has an abundance gradient (Watanabe et al. 1997; Fujita 
et al. 1996; Tamura et al. 1996; Finoguonov et al. 1999). It 
is possible that mergers that could be responsible for the 
disruption of the cooling flow in clusters can also destroy 
abundance gradients. The degree to which the abundance 
gradient remains intact may be an indicator of how recent 
and/or violent the merger was. All four non-cooling flow 
clusters are merger candidates. Of the five cooling fiow 
clusters with significant gradients, only A85 is a merger 
candidate. Of the three cooling fiow clusters with only 
weak evidence for an abundance gradients two are recent 
merger candidates. Oegerle & Hill (1994) conclude that 
the velocity distribution of galaxies within A1795 suggest 
that the cluster underwent a merger with a relatively small 
subcluster such that the cooling flow was not disrupted. 
Oegerle, Hill, & Fitchett (1995) also find evidence for sub- 
clustering from the galaxy velocity distribution in A2142. 
ROSATPSPC and HRI images of A2142 also suggest that 
a merger is in process (Henry & Briel 1996). It is possi- 
ble that mergers have mixed the gas of these two clusters 
enough to decrease the significance of an abundance gra- 



dient without destroying the cooling fiow. A2029 does not 
show evidence for a merger, although the lack of a signifi- 
cant abundance gradient may be the result of poor photon 
statistics rather than a lack of a gradient. Further studies 
of A2029 will be needed to confirm this. 

As suggested by Fabian et al. (1994) and Allen & Fabian 
(1998), the presence of abundance gradients in cooling 
flow clusters provides a reasonable explanation as to why 
cooling fiow clusters have higher emission-weighted metal 
abundances than their non-cooling fiow counterparts. If 
the cooling gas at the center of clusters has a higher 
abundance than the outer regions, the strongly peaked X- 
ray surface brightness distribution characteristic of cooling 
flow clusters will preferentially weight the global emission- 
weighted abundance towards the central value, with the 
low-abundance outer regions contributing relatively little 
to the global value. Such an effect is not present in non- 
cooling fiow clusters owing to the broader, less-peaked X- 
ray surface brightness distribution of these clusters. But 
is this the only reason for the higher measured metallicity 
in cooling fiow clusters, or is there really an abundance 
difference between cooling flow and non-cooling flow clus- 
ters? To test this, we determined abundances in two re- 
gions for each cluster, a circle with a radius of 0.075ryiriai 
and an annulus with inner and outer radii of 0.075ryiriai 
and 0.173ryiriai- The outer radius was set by the need 
to extend the annulus no further than 9' for any cluster, 
and the inner radius was set to be larger than the cooling 
radius for any of the cooling flow clusters. 

The results are shown in Table [|. Inside of 0.075rviriai 
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Fig. 2. — Normalized abundance profiles for all 12 clusters in the sample versus radius in units of the virial radius. Here, the error bars 
represent the la errors. A significant radial gradient in the abundance is evident. The solid line represents the composite abundance profile 
of 18 realizations of the cluster simulations of Metzlcr & Evrard (1997). 



the cooling flow subsample had an average abundance of 
0.42 ± 0.06, while the non-cooling flow subsample had an 
average of 0.33 ± 0.04. In the outer region, the cooling 
flow subsample had an average abundance of 0.30 ± 0.02, 
while the non-cooling flow subsample had an average of 
0.24±0.03. Thus, there appears to be an abundance gradi- 
ent in non-cooling flow clusters as well as cooling flow clus- 
ters. The presence of abundance gradients in non-cooling 
flow clusters (especially in hot clusters) would seem to rule 
out claims that previously observed abundance gradients 
were artifacts of incorrect spectral modeling due to either 
the presence of cooling gas at the center of clusters or un- 
certainties in modeling the Fe-L line complex for cooler 
clusters. 

Although the average abundances of the cooling flow 
and non-cooling flow clusters are marginally consistent 
within the errors within the cooling region, the data sug- 
gests that there is a discrepancy between the abundance of 
cooling and non-cooling flow clusters outside of the cooling 
radius. Of course, any conclusions drawn from a subsam- 
ple of four clusters should be viewed with caution. A larger 
sample of clusters will be necessary to confirm if the abun- 
dance differences between cooling flow and non-cooling 
fiow clusters are real. This result is in the opposite sense 
expected if higher mctallicity gas from the cooling flow re- 
gion has been mixed in with the lower metallicity gas of 
the outer regions as the result of a merger in a present- 
day non-cooling flow cluster. This trend also appeared in 



a sample of clusters observed by ASCA. Fukazawa et al. 
(1998) analyzed 26 clusters with temperatures of at least 
3.0 keV, and found global abundances after having excised 
the inner cooling regions. The average abundance for the 
16 cooling flow clusters was 0.29 ± 0.07, while the average 
for the ten non-cooling flow clusters was 0.23 ± 0.05. Al- 
though at a lower significance than our BeppoSAX sample, 
the trend is in the same direction. Fukazawa et al. (1998) 
did not claim to extract their spectra in units of the virial 
radius, which is the unit of distance that allows the most 
direct comparison of radial properties for clusters of vary- 
ing sizes and temperatures. This might wash out some of 
the information on the abundance gradient when the pro- 
files are combined, and lead to a result of lower statistical 
significance. 

Allen & Fabian (1998) present the interesting scenario 
in which metallicity gradients can be explained (at least 
in part) by a significant fraction of the metals in clusters 
residing on grains in cluster cores. Since the lifetime of 
a grain to sputtering in the hot ICM is inversely propor- 
tional to the electron density of the hot gas (e.g., Draine & 
Salpeter 1979), more metals would be released into the gas 
phase in the high density cooling flow regions of clusters 
than in lower density regions, providing the dust grains 
are > 10 microns. By extending this argument outside 
the core region, it might be possible to explain the differ- 
ence in metallicity at large radii between cooling flow and 
non-cooling flow clusters. If the density of gas outside the 
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-0.03 
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n 4'^+005 
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-50+003 
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A2256 


7-lt^1 
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9.9ir7 
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t. -1+0.6 


2A0335+096* 


q 9 + 0.08 

•J-^-o.os 



0.37^ 



0.35 
0.28 
0.28 
0.24 
0.40 
0.41 



-0.06 
+0.02 
-0.02 
-0.04 
-0.04 
-0.05 
-0.05 
-0.05 
-0.05 
-0.08 
-0.08 
-0.03 
-0.03 



1.17/282 
1.20/259 
0.97/211 
0.98/141 
1.02/292 
1.05/440 
1.03/419 
1.06/238 
1.04/273 
0.93/265 
0.95/165 
1.13/257 



0.45^ 



0.36^ 



0.22 

0.23^ 

0.30^ 



* cooling fio-w cluster 

"Best-fit global temperature of BeppoSAX data from Irwin 



Brcgman (2000), except for A3562 which is presented here for the first time. 



core is greater in cooling flo-w clusters than in non-cooling 
flo-w clusters, more metals -will be released into the gas 
phase for the former. A difference in density among cool- 
ing and non-cooling flo-w clusters outside of 200 kpc -was 
observed by White et al. (1997), -who deprojected the sur- 
face brightness profiles of 207 clusters observed -with Ein- 
stein Observatory and found that the electron densities of 
clusters "with cooling rates greater than 10 Mq yr~^ -were 
larger than for clusters -with cooling rates less than 10 Mq 
yr^^ out to a radius of 1 Mpc. For a 7 keV cluster, our 
0.075 — O.nSryiriai bin corresponds to 240-560 kpc. At 
a radius of 450 kpc. White et al. (1997) found that the 
cooling flo-w clusters had an electron density of ~7 x 10^^ 
cm^3^ -whereas the non-cooling flo-w clusters had an elec- 
tron density of ^3 x 10^'' cm^^^ por a grain lifetime of 
2 X 10^ a/rie yrs (-where a is the size of the grain in microns 
and Tie is the electron density of the hot gas in cm^^)^ 
grains -with a size of around 2 microns -would survive in 
non-cooling flo-w clusters for a Hubble time while being 
sputtered in cooling flo-w clusters. While rather large, dust 
grains of this size are not unreasonable. In addition, this 
argument assumes that the dust -was deposited into the 
ICM at the time the cluster -was created. In reality, dust 
is deposited into the ICM throughout the history of the 
cluster, allo-wing smaller dust grains to survive in lo-w den- 
sity non-cooling flo-w clusters -while still being sputtered in 
high density coling flo-w clusters. This in turn can lead 
to a higher abundance in the outer regions of cooling flo-w 
clusters since more of the grains have been sputtered away 
owing to the higher electron density. 

However, more recent determinations of the density pro- 
files of clusters determined with ROSAT PSPC data in- 
dicate that the non-cooling fiow clusters in our sample 
do not have a lower electron density in the 0.075-0.173 
Tviriai region than the cooling flow clusters (J. Mohr, pri- 
vate communication). At 10% of the virial radius, the 
average electron density for the non-cooling flow clusters 
A2256, A2319, and A3266 was 1.80 ± 0.20 x 10"^ cm'^, 



while for the cooling flow clusters A85, A496, A1795, 
A2029, A2142, A2199, and A3562 the average electron was 
1.92 ± 0.46 X 10^3 cm^3 xhus, it does not appear likely 
that the difference in abundances outside the cooling ra- 
dius can be the result of an increased sputtering of grains 
in cooling flow clusters. 

Since non-cooling flow clusters are usually merger candi- 
dates, one possible explanation for the difference in abun- 
dances between cooling and non-cooling flow clusters is 
the mixing of gas via mergers from larger (> QAlbrmriai) 
radii rather than from gas interior to this region. If the iron 
abundance of the gas continues to decline with increasing 
radius beyond Q.2rmriai as the simulations of Metzler & 
Evrard (1997) imply, the mean metallicity in our 0.075- 
0.173 ryiriai bin will be lowered upon mixing with this low 
metallicity gas during a merger. A2142, a cooling flow 
cluster that is thought to be undergoing a merger and 
shows only weak evidence for an abundance gradient, also 
has the lowest metallicity in the 0.075-0.173 ryiriai bin 
among cooling flow clusters in our sample, consistent with 
the merger-mixing theory. However, A1795, which is also 
believed to have recently experienced a merger has a high 
metallicity value in the 0.075-0.173 ryiriai bin, although 
the uncertainty is rather large. A larger sample of clus- 
ters with abundances determined out to larger radii will 
be needed to explore this topic further. 

4. CONCLUSIONS 

We have analyzed 12 clusters with the MECS instru- 
ment onboard BeppoSAX and flnd evidence for a nega- 
tive abundance gradient in most of them. This gradient 
is present in all eight cooling flow clusters, and to a lesser 
significance in the non-cooling flow clusters. Since kT > 5 
keV for nine of the 12 clusters, this work extends the ubiq- 
uity of abundance gradients to hot clusters. The slope of 
the abundance gradient is in good agreement with cluster 
simulations performed by Metzler & Evrard (1997). In our 
sample cooling flow clusters had a higher metallicity than 



non-cooling flow clusters both in the inner regions of the 
cluster and the outer regions out to several hundred kpc; 
further work with a larger sample will be needed to confirm 
this result for clusters in general. Outside the cooling re- 
gion, it seems unlikely that this difference can be explained 
by the more efhcient sputtering of dust grains in cooling 
flow clusters, since the densities of cooling flow and non- 
cooling flow clusters are similar in this region. The mixing 
of low metallicity gas during a merger may be responsible 
for this difference. 

We thank the anonymous referee for many useful sug- 
gestions and comments. We thank Chris Metzler for mak- 



ing his theoretical abundance profiles available to us. We 
also thank Joe Mohr for kindly providing us his ROSAT 
PSPC density profiles. JAI thanks Renato Dupke and 
Gus Evrard for many useful comments and conversations. 
This research has made use of data obtained through the 
BeppoSAX Science Data Center and the High Energy As- 
trophysics Science Archive Research Center Online Ser- 
vice, provided by the NASA/Goddard Space Flight Cen- 
ter. This work has been supported by Chandra Fellowship 
grant PF9-10009, awarded through the Chandra Science 
Center. The Chandra Science Center is operated by the 
Smithsonian Astrophysical Observatory for NASA under 
contract NAS8-39073. 



REFERENCES 



Allen, S. W., & Fabian, A. C. 1998, MNRAS, 297, L63 

Anders, E., & Grevesse, N. 1989, Geochimica et Cosmochimica Acta, 

53, 197 
Boella, G., Chiappetti, L., Conti, G., Cusumano, G., Del Sordo, S., 

La Rosa, G., Maccarone, M. C., Mineo, T., Molendi, S., Re, S., 

Sacco, B., & Tripiciano, M. 1997, A&A, 122, 327 
D'Acri, F., De Grandi, S., & Molendi, S. 1998, in The Active X-ray 

Sky: Results from BeppoSAX and RXTE, ed. L. Scarsi, H. Bradt, 

P. Ginrnmi gnri T7 F^nrt^ ( A matorrlgm ■ Plgmncry FjRI 



Kaastra, J. S., Bleeker, J. A. M., & Mewe, R. 1998, in The Active 
X-ray Sky: Results from BeppoSAX and RXTE, ed. L. Scarsi, H. 
Bradt, P. Giomini, and F. Fiore (Amsterdam: Elsevier), 567 

Kikuchi, K., Furusho, T., Ezawa, .H , YflmflSRki N., Ghashi, T., 

Fukazawa, Y., & Ikebe, Y. 1999, astaxtpii/aaOilil . 

Metzler, C. A., & Evrard, A. E. 199r"^ ' 



^stro-pli/'9 1 1U324 



Metzler, C. A., & Evrard, A. E. 1994, ' ApJ, 43V, 564 

Mohr, J. J., Mathiesen, B., & Evrard, A. E. IQQQ ApT .^17 fi?? 



De GrAndi, S., & Molendi, S. 1999a, A&A, 351, L45 
De Gr andi, S., fc Molendi, S. 1999b, ApJ, 52V, L25 

nraino VI. T Rr Bglpofor P, P. 1 Q7Q ApT O'iA 77 



A/fnlonrH B 1 QQ« Rppp^ ?/) V Tor'tinirgl Report, ftp://www.sdc. asi.lt/pub/sax/d 



Dupke R. A., & White, R. E., III. 2000a, ApJ, 528, 139 

Dupke, R. A., fc White, R. E., III. 2a0ab, ApJ, in press 

Ettori, S., & Fabian A. C. 1999, MNRAS, 305, 834 

Evrard, A. E., Metzler, C. A., & Navarro, J. F. 1996, ApJ, 469, 494 

Ezawa, H., Fukazawa, Y., Makishima, K., Ohashi, T., Takahara, F., 

Xu, H., & Yamasaki, N. Y. 1997, ApJ, 490, L33 
Fabian, A. C., Crawford, C. S., Edge, A. C, & Mushotzky R. F. 

1994, MNRAS, 267, 779 



Finogupnov, A., David, L. P., & Ponman T. J. 1999, astro- 
ph/! l908150| 

Fujita, v., Koyama, K., Tsuru, T., & Matsumoto, H. 1996, PASJ 
48, 191 

Fukazawa, Y., Makishima, K., Tamura, T., Ezawa, H., Xu, H., Ikebe 
Y., Kikuchi, K., & Ohashi, T. 1998, PASJ, 50. 187 

Gaetz, T. J., Salpeter, E. E., Shaviv, G 1987, ApJ, 316, 530 

Gunn, J. E., & Gott, J. R. 1972, ApJ, 176, 1 

Henry J. P., & Briel, U. G. 1996, ApJ, 472, 137 

Ikebe, Y., et al. 1997, ApJ, 481, 660 

Irwin, J. A., & Bregman, J. N. 2000, ApJ, in press 

Ishimaru, Y., & Arimoto, N. 1997, PASJ, 49, 1 



extended_sources.ps. gz 
Molendi, S., fc De Grandi , S. 1999, A&A, 351, L41 
A/TnlianrH S Dr. Grandi, S., Fusco-Femiano, R. 



ph/0d03060| 
)lcndi, S., De 



2000, astro- 



Ivlolcndi, S., De Grandi, S., Fusco-Femiano, R. , Colafrancesco, S., 

Fiore, F., Nesci, R., & TambureUi, F. 1999, ApJ, 525, L73 
Oegerle, W. R., & Hill, J. M. 1994, AJ, 107, 857 
Oegerle, W. R., Hill, J. M., & Fitchett M. J., AJ, 110, 32 
Ostriker, A., Yahil, J. P. 1973, ApJ, 185, 787 
Peres, C. B., Fabian, A. C, Edge, A. C, Allen, S. W., Johnstone, R. 

M., & White, D. A. 1998, MNRAS, 298, 416 
Sarazin, C. L., Wise, M. W., & Markevitch, M. L. 1998, ApJ, 498, 

606 
Takahashi, T., Markevitch, M., Fukazawa, Y., Ikebe, Y., Ishisaki, 

Y., Kikuchi, K., Makishima, K., & Tawara, Y. 1995, ASCA 

Newsletter, #3 (NASA/GSFC) 
Tamura, T., et al. 1996, PASJ, 48, 671 
Watanabe, M., Yamashita, K., Kunieda, H., & Tawara, Y. 1997, in 

X-ray Imaging and Spectroscopy of Cosmic Hot Plasmas, ed. F. 

Makino & K. Mitsuda (Tokyo:IJnii,'ersal Academy), 131 
White, D. A. 1999, a.stro-ph/9909467 
White, D. A., Jones, C, Forman, W. 1997, MNRAS, 292, 419 



